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Abstract 

In this paper we prove a strong maximum principle for certain parabolic systems of equations. In 
particular, our methods place no restriction on the regularity of the boundary of the convex set in which 
the system takes its values, and therefore our results hold for any convex set. We achieve this through 
the use of viscosity solutions and their corresponding strong maximum principle. 

Introduction 

In this paper we prove a strong maximum principle for certain parabolic systems. To be precise, let X C 1" 
be open and connected, and set 17 = Ix (0,oo). Suppose u G C 2,1 (fl;M. k )nC(Cl;'R k ) satisfies a parabolic 
system of equations of the form 

— = D(x, t, u) 2^ (Hj{x, t) q—q- + ^ Mi(x, t, u)— + <f>(x, t, u). (1) 

i,j % 3 i 

where the aij are real-valued, <p takes values in R fc , and D(x, t, z) and each of the Mj(x, t, z) take values in 
the space of k x k matrices. 

We make the following regularity assumptions: As a function of z, we assume that 4>{x, t, z), D(x, t, z) 
and each of the Mi(x,t, z) are Lipschitz continuous uniformly for (x,t) in compact subsets of 17. We 
assume also that each of the aij(x, t) is locally bounded in Q and that D(x, t, z) and each of the Mj(x, t, z) 
are locally bounded in 17 x M. k . Finally, we assume that the matrix {aij(x,t)} is symmetric and locally 
uniformly positive definite in 17 and the matrix D(x, t, z) is locally uniformly positive definite in 17 x R fc . 

Next, suppose that K is a closed, convex subset of M fe which is compatible with <p, D, and the M% in the 
following sense: For all (x, t) G 17, v G dK and each vector v which is inward pointing at v (See § 1 for 
the definition of "inward pointing at v."), we have that cj)(x, t, v) ■ v > and that v is a left eigenvector of 
D(x, t, v) and each of the Mj(x, t, v). 

We will show that if 

u(x, t) e K for all (x, t) G 17, (2) 

then 

u(xq, to) ^ 9K for some (xq, t ) G 17 => u(x, t) G dK for all (x, t) G X x [0, to], (3) 
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and we will refer to the implication in (O as the strong maximum principle for solutions to (Q]). Under 
circumstances when solutions to (GQ) satisfy a weak maximum principle of the form 

u(x, t) G K for (x, t) edn => u(x, t) G K for (x, t) G 0, (4) 

it is obvious that © can be replaced by 

u(x, t) G K for (x, t) G 90. 

In his 1975 paper [6], H. Weinberger considers the case where D = I and the Afj are real-valued 
functions. In this case, Weinberger states and proves the strong maximum principle for parabolic systems 
under a regularity condition on the convex set K which he calls the "slab condition". In their 1977 paper HI, 
K. Chueh, C. Conley, and J. Smoller prove, under mild regularity assumptions, a weak maximum principle 
of the form ((U) for the system of partial differential equations (O, under the same compatibility condition 
we have assumed on K, <fi, D, and the Mj. In his 1990 paper [5], X. Wang, gives a geometric proof of the 
strong maximum principle for dD, following Weinberger's arguments, in the case that the boundary of K is 
C 2 . In this case, the distance function, d, to the boundary of K is C 2 in K (at least near the boundary), in 
which case he proves the strong maximum principle by applying the classical strong maximum principle to 
d(u(x, t)). 

The argument presented in this paper is similar to X. Wang's; we too will apply the classical strong 
maximum principle to d(u(x,t)). However, we are able to remove all regularity requirements on dK by 
showing that even if d(u(x, t)) is not twice differentiable, it is still a viscosity super solution of the appro- 
priate partial differential equation (See [2] for an introduction to the theory of viscosity solutions). We can 
then apply a strong maximum principle for viscosity solutions to parabolic partial differential equations to 
complete the proof of ([3]). The strong maximum principle we need is provided by F. Da Lio in |3], wherein 
Da Lio extends the classic argument given by L. Nirenberg in [4] to the case of viscosity solutions. 

1 Proof of© 

We begin with the following definitions: 

Definition 1.1. Given a convex set K and a boundary point v G dK, a function £ : K — > R is called a 

supporting affine functional of K at v if £{z) > for all z G K, £{v) = 0, and |V-£(z)| = \V£\ = 1 (That 
is, £(z) is the distance to a supporting hyperplane of K at v, which, by abusing notation, we also denote by 
£). 

We say that a function £ : K — > E is a supporting affine functional of K if it is a supporting affine 
functional of K at some v G dK. 

Definition 1.2. For v G dK, a vector v is an inward pointing vector at v G dK if there exists a supporting 
affine functional £{z) at v such that V£ = v (That is, if there exists a supporting hyperplane £ of K at v 
whose "inward pointing" unit normal vector is v). 

Note that in our definitions we include the assumption that an inward pointing vector has unit length and 
that a supporting affine functional has unit-length gradient. For geometric reasons, these are nice assump- 
tions to make as they allow for the intuitively appealing interpretation of £{z) as the distance to a supporting 
hyperplane £. 

Given v G dK, we use L v to denote the set of all the supporting affine functionals £ of K at v. By 
well known results (e.g., the Hahn-Banach Theorem), L v ^ for every v G dK. For z G K we let 
d(z) = ini{\z — v\ : v G dK} denote the distance from z to the boundary of K. 

First we present a lemma from convex analysis. 
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Lemma 1.3. For any convex set K, 

d(z) =inf t(z), 

where the infimum is taken over all supporting affine junctionals £. Furthermore, for each z G K, there is a 
(possibly non-unique) point v G dK and an £ G L v such that 

d(z) = \z — v\= £(z). 

In fact, L v consists of just this one supporting affine functional £. 

Proof. It is easy to see geometrically that d(z) < £(z) for each £: Given £, l(z) = \z — z\, where z is the 
projection of z onto the supporting hyperplane determined by I. The line from z to z intersects dK at some 
point w. Thus 

d(z) < \z — w\ < \z — z\ = £(z). 

Next we show that the infimum over I is in fact achieved. Using B(z, r) to denote the closed ball with 
radius r centered at z, we have that d(z) = ini{\z — w\ : w G dK} = inf{|z — w\ : w G 8KnB(z, 2d(z))}. 
By compactness, this infimum is achieved and we have d(z) = \z — v\ for some v G dK. 

Take £ G L v . We have shown that d(z) < £(z). Since v G £, we also have that £[z) = dist(z,£) < 
\z — v\= d(z). And so it follows that d(z) = £(z). 

So for all £ G L v we have that £(z) = \z — i;| which means that the line from z to v is normal to the 
hyperplane £. As there is only one hyperplane with this normal which touches v, we have that L v is just the 
singleton {£}. 

□ 

In the remainder of the paper we will use the notation d(x, t) for d(u(x, t)) and £(x, t) for £(u(x, t)). 
The following is the key result proved in this section. 

Theorem 1.4. d(x, t) satisfies, in the viscosity sense, a parabolic equation of the form 
dd d 2 d sr-^ n dd - 

et -^ a ^dx~dx- -L^ f t+^ f ) d ^ o. 

where j(x, t) > in Q, each of the ctij,{3i, and 7 are locally bounded in £1, and {otij} is locally uniformly 
positive definite in £1. 

Proof. Call a quadruple (x, t, v, t) nice if 

1. (x, t) G fl,v G dK, and I is a supporting affine functional of K at v. 

2. d(u(x, t)) = \u(x, t) — v\. 

By the previous lemma, we know that for each (x, t) we can find at least one v G dK satisfying the 
second condition and this in turn will determine a unique choice of £ satisfying the first condition. We then 
also have that £(u(x, t)) = |u(x, t) — v\. 

For any nice quadruple we have that 
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V£-u t 

V£- [D{x,t,u) s ^a ij {x,t)\x XiX] Mi(x,t, u)u Xi + 4>(x,t,u)} 

V£- [D(x,t,u)^2aij(x,t)u XiXj Mi(x,t, u)u Xi + (f)(x,t,u)} 

-VI- [D(x,t,v)^2aij(x,t)u XtXj +^2Mi(x,t,v)u Xi +4>(x,t,v)} 

i,j i 

+V£- [D(x,t,v)^2aij(x,t)u XiXj + ^ M { (x, t, v)u x . + <j)(x,t,v)} 



> 



u — v\ 



c(x, t) || ^ ajjjx, t)ua, <3;i . 1| + y^mj(x,t) \\u Xi \\ + p(x,t) 

t, v, VI) aij(x, t)(V£- u XiXj ) + Y, x i( x > *> v > • + 



c(x, t) || ^ ajjjx, t)^^ || + y^mj(g;,t) ||u x J| +p(x,t) 
»J » 

+yu(a;, t, v, VI) ay (a:, + A* (a, t, u, V^)4 4 , 

where c(x,t), nii(x,t), and p(x,t) are the Lipschitz constants for D(x,t,-), Mi(x,t, •), and (j)(x,t,-), 
respectively, and ^ and Aj are the eigenvalues for D and Mi at i> corresponding to the left eigenvector 
VI (recall that VI is inward pointing at v). Note that in the inequality step we have used the fact that 
4>(x, t, v) ■ Vi > which follows from the fact that VI is inward pointing at v and from our compatibility 
assumption on <j) and K. 

Next, let 7 (x, t) = c(x,t)\\Y^ij aij(x,t)u XlXj \\ +J2i m i( x i t ) \\ u Xi\\ + p(x,t) ■ It follows from our 
regularity assumptions on D, Mi, and (j> that c(x,t), rrii(x,t), and p(x,t) are locally bounded. Since 
u G C 2,1 (Q), its first and second spatial derivatives are locally bounded, and so, since the ay(z;, i) are 
locally bounded, it follows that j(x,t) is locally bounded. As each of the Lipschitz constants is positive, it 
is clear that 7(2;, t) > 0. 

Thus we have shown that for each nice quadruple (x, t, v, £), 



/i(2?,t,w,V^) J^ayfot)^ -^A;(x,M,W)4 +-y(x,t)£> 0. 



(5) 



As remarked at the outset of this proof, given (x,t) G SI we can always find awe dK and £ to form a 
nice quadruple (x, t, v, £). In general there will be more than one way to extend (x, t) into a nice quadruple, 
but for our purposes we only care that an extension is possible. 

To avoid the axiom of choice, we now describe a method of choosing an extension: Given (x,t) G Q we 
need that v G dK be such that \u(x, t) — v\ = d(u(x, t)). The set of v which satisfy this relation is a closed 
and bounded and hence compact subset of M fc . We first look at v in this set with smallest first component. If 
there is a unique such v we choose it. Otherwise, among those v we look for the one with smallest second 
component and we continue this algorithm until we find a unique v. We denote this v by v^ x ^ to make clear 
its dependence on (x, t). Lastly, once we have v, £ is uniquely determined (see Lemma [T3T ) and we denote 
it by 



4 



Next we let j~i(x, t) = /x(x, t, v/ Xtt -\,V£/ X) t)) and Xi(x, t) = \i(x, t, vi x>t \,V£i x ^))- We claim that jl(x, t) 
and the Aj(x, t) are locally bounded. This is true as for any compact set C C Vl, 

sup £t(x,t) < sup sup n(x,t,v,V£) 

(x,t)eC (x,t)eC {v:\u- v\=d(u)} 

< sup sup < oo, 

(x,i)eC {v.\u-v\=d(u)} 

where \\D\\ denotes the operator norm of D which is locally bounded by our regularity assumption on D. 
The same argument works for the Aj. 

We next claim that fi(x, t) is uniformly bounded away from on compact sets. This is true as for any 
compact set C C ft, 

inf u(x,t) > inf inf u(x,t,v,'V£) 

(x,t)ec (x,t)ec {v.\u-v\=d(u)} 

> inf inf Ai(D(x,t,u)) > 0, 

(x,t)eC {v:\u-v\=d(u)} 

where we denote by A\(D) the smallest eigenvalue of the positive definite matrix D. Here we have used 
our assumption that D(x, t, v) is uniformly positive definite on compact sets. 
Finally we reach the crux of our argument. We claim that d solves 

d t - fL(x,t) J ^2a ij (x,t)d XtXj - } ]Xj(x,t)d x . + j(x,t) > 

i,3 * 

in the viscosity sense. 

Suppose that ift G C°°(fl) touches d from below at the point (x,t), i.e. 

ip(x, i) = d(x, t), and ip(x, t) < d(x, t) in a neighborhood of (x, i). 

We extend the point (x, i) to the nice quadruple (x, i, v, £ ft,£/ x £))• Since by Lemma [T31 d(z) < £,~ ^(z), 
we then have that 

tp(x,t) = £,~fa(x,i), and ip(x,t) < £, x ^{x,t) in a neighborhood of (x,i). 
Therefore, ip also touches the function £,* ^ from below at (x,i), and so 

d d - - 

of = ^(4,f)(*.*). and A ^ < AJ {Xjt) {x,t). 

It then follows from (0) (and recalling the definitions of jj, and Aj) that 

Since this is true for all points (x, t) and all smooth functions tp touching d from below at (x,t), we have 
shown that d(x, t) solves 

dt - fi(x, t) 

t)d XiX . - ^2 + 7(s> t) d > 

i,j i 

in the viscosity sense. 

The theorem is now proved by letting = A, and ctij = jldij. Note that {a^} is uniformly positive 
definite on compact sets as } is uniformly positive definite on compact sets and fl is uniformly bounded 
away from on compact sets. 

□ 
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In order to complete the proof of (O, we need the following version of the strong maximum principle 
for supersolutions in the viscosity sense: 

Theorem 1.5. Suppose that f E C(J); [0, oo)) satisfies 

% - £ ««(«. ^ " £ «<*• + 7(I - f)/ - 

/« f/ie viscosity sense. Suppose further that in Q., 7 > 0, eac/j of the au, and 7 are locally bounded, and 
{ctij} is locally uniformly positive definite. Then 

/(xq, to) = 0/or some (xo, to) £ ^ ==> f( x , t) = Ofar all (x, t) £ Cl with t < to. 

Note that by applying this theorem to / = d > and using Theorem 11.41 we get © as an immediate 
consequence. 

Proof. By linearity, we have that — / solves 

% -£MM)^- "EftM-g +7(M)/ < (6) 

in the viscosity sense, and by our non-negativity assumption on /, we have that — / achieves a non-negative 
maximum at (xo, to)- If — / £ C 2,1 (S7; [0, 00)) were a classical solution to ©, it would immediately follow 
from the classical strong maximum principle that — / = / = for all (x, t) E with t < to- 

As — / only solves © in the viscosity sense, we instead rely on the extension of the proof of the classical 
strong maximum principle given by Nirenberg in [4] to the case of viscosity solutions, provided by F. Da 
Lio in 0. Da Lio's assumptions are more general than we need and apply to viscosity solutions of more 
general PDE of the form 

G(x,tJJ t ,D x f,D 2 xt f)<0 

under the assumption that G is continuous. We don't have continuity for the coefficients of our PDE ©, 
and so our PDE does not directly fall under the assumptions of Da Lio's result. Nevertheless, because our 
PDE © is of such a simple form, Da Lio's method of extending Nirenberg 's proof can still be applied to 
our PDE, and we achieve the same result that — / = / = for all (x, t) E ft with t < to- □ 



2 Acknowledgments 

I thank my advisor, Professor Daniel W. Stroock, for his numerous helpful conversations and suggestions. I 
also thank Professor David Jerison for his helpful advice. 



References 

[1] K. Chueh, C. Conley, and J. Smoller, Positively invariant regions for systems of nonlinear diffusion 
equations, Indiana University of Mathematics Journal 26, 373-392 (1977). 

[2] M. Crandall, H. Ishii, P. L. Lions, User's guide to viscosity solutions of second order partial differential 
equations. Bull. Amer. Math. Soc. 27, 1-67 (1992). 

[3] F. Da Lio, Remarks on the strong maximum principle for viscosity solutions to fully nonlinear parabolic 
equations, Commun. Pure Appl. Anal. 3 (3) 395-415 (2004). 



6 



[4] L. Nirenberg, A strong maximum principle for parabolic equations, Comm. Pure Appl. Math., 6, 167- 
177 (1953). 

[5] X. Wang, A remark on strong maximum principle for parabolic and elliptic systems, Proceedings of the 
American Mathematical Society, Vol. 109, No. 2, pp. 343-348 (Jun.,1990). 

[6] H. Weinberger, Invariant sets for weakly coupled parabolic and elliptic systems, Rend. Mat. (7) 8, 295- 
310(1975). 



7 



